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Quasinormal mode (QNM) of black holes has been an intriguing subject of discussions for the last three 
decades Q, O, [|. It is believed that QNM is a characteristic sound of black holes which could lead to the 
direct identification of the black hole existence through gravitational wave observations to be realized in the 
near future Q, In addition to its potential astrophysical interest, theoretically QNM was believed as a tool 
to learn more about black hole and was even argued as a testing ground for fundamental physics. It was found 
that the QNMs of anti-de SitterfAdS) black holes have direct interpretation in terms of the dual conformal 
field theory (CFT) 3, U, la, lO, H 13, l2( • This could serve as a support of the AdS/CFT correspondence. Attempts 
of using QNMs to investigate the dS/CFT correspondence have also been proposed [lOj. Recently it was argued 
that QNMs might reflect the possible connection between the classical vibrations of a black hole spacetime 
and various quantum aspects by relating the real part of the QNM frequencies to the Barbero-Immirzi(BI) 
parameter, which was introduced by hand in order that loop quantum gravity reproduces correctly the black 
hole entropy [11 



12]. But tin 1 direct connection lias not been found in AdS black hole background Q. 



It is of great interest to investigate whether QNM can reflect more physics of black holes. Recently some 



14 



15, 



16 1 . 



indications have been found that black hole phase transitions can show up in the QNM spectrum 
This is interesting since it might be the first phenomenon telling us the existence of the phase transition 
in black hole physics. In order to examine whether the QNM is an effective probe of phase transitions, we 
need to investigate more general black hole configurations with more general field perturbations. Jing et al 
171 ] computed the QNM of Reissner-Nordstrom (RN) black hole and claimed that they found the second 



19(, it was argued that the result in 



order phase transition point predicted by Davies 18| where the heat capacity appears singular. Bowever in 



171 ] might probably be a numerical coincidence and the conjectured 



correspondence between QNM and Davis's phase transition does not straightforwardly generalize to Kerr or 
Schwarzschild AdS metrics. Whereas calculations in li| cannot rule out the relation between dynamical and 
thcrmodynamical properties of black holes but suggest that such a relation is non-trivial. 

In the study of black hole phase transition, there have been a lot of debate on Davies' phase transition 
point. Around Davies' point, the black holes' event horizons do not lose their regularities and internal states 
of black holes do not change significantly. It was considered more reasonable that the phase transition of 
black hole occurs when a nonextremal black hole approaches its extremal limit, since all second moments of 

n 

non-equilibrium fluctuations diverge there characterizing the second order phase transition [20|. Moreover 
extremal black holes are very different from nonextreme holes. Extremal holes just have superradiation but 
without Bawking radiation since its Bawking temperature vanishes. The extremal black holes' geometric 



structures are also different from their nonextremal counterparts, since the singularity will be naked beyond 
the extremal limit. Davies' points for the divergence of the heat capacity were later proved in fact as turning 
joints which are related to the changes of thermodynamical stability, especially in the canonical ensemble 



21 1 . Whether the thermodynamical stability can be reflected from the dynamical stability such as the QNM 



spectra has been under heated discussions recently (for a review see [22|). For a black brane solution, it 



23] . But whether this 



was argued that the thermodynamical stability is often related to dynamical stability 
correspondence is profound in all spacetimes is still not clear. 

In this work we are going to study the QNM in the background of charged Kaluza-Klein(KK) black hole 
with squashed horizons [2j| and investigate whether the QNM spectra can reflect the Davies' points of ther- 



modynamical stability in this background. In 25j, the thermodynamic properties of the KK black hole are 
discussed and compared to its undeformed five-dimensional RN black hole counterpart. The Davies' ther- 
modynamical stability point was found where the heat capacity diverges when the horizon of the black hole 
crosses the critical value. We will examine the scalar perturbation in this background and see whether the 
thermodynamic stability can be reflected in the QNM spectrum. 

The five-dimensional charged KK black hole with squashed horizons is described byj^J] 

1-2 1 \ 2 2 

ds 2 = -f{r)dt 2 + J^j- dr2 + ^jk{r)dSf + T — {d^j + cos Bd^f , (1) 
where dQ 2 = dO 2 + sin 2 9dip 2 is the metric of the unit sphere and 



(r 2 - r\){r 2 - r 2 _) ^ (r*, - r 2 ,)^ - r 2 _) 

{rlo - r 2 ) 2 



/W - ^3 -> Kr) = V °° ,, 2 + - • (2) 



Here < 6 < n, < (f> < 2ir and < ip < An are Euler angles. The gauge potential is given by 



V3 



r + r- 



As in the RN black hole, the coordinate singularities r — r + and r — r_ correspond to the outer and inner 
horizons of the black hole, respectively. r w is the spatial infinity. In the parameter space < r_ < r + < r^, 
r is restricted within the range < r < Too. The shape of black hole horizon is deformed by the parameter 
k(r+). 

In the metric (1), the intrinsic singularity is just the one at r = 0. This can be seen by introducing a new 
radial coordinate p as 



with 



2 _ fco 2 



fco = k(r = 0) 



(•5) 



At the spatial infinity r — > roo,P — > oo. Thus in the new coordinate, p varies from to oo when r changes 
from to Too. The metric ([T]) can be rewritten as 



ds 2 = -F{p)dr 2 + ^j$dp 2 + p 2 K 2 (p)dn 2 
F(p) 



4K 2 (p) 



(dip + cos 6d(f>) 2 , 



with 



Po 



K 2 (p) = l + ^ 



(6) 



(7) 



P J \ P J P 
Here we have defined the proper time r = 2pot/r cxl for the observer at infinity. The mass and charge for this 
squashed KK black hole are defined by 



V37T7-0O 

<2 = g y/P+P-' 



(8) 



where p± = Po r ±/( r ^o — r ±) are the outer and inner horizons of the black hole in the new coordinate 



The Hawking temperature and entropy of the black hole can be expressed as 

— : ~2 J2 „2 



25] 



T, 



H 



P+ ~ P- 

4np 2 + y p + + p *„, + loo 



9 2 



S = 4tt 2 (p + ) ^ (po + p _) 2 + p _) = 



7T 2 ri r 2 — r 2 . 



2 - r 

Several limits of this charged squashed KK black hole were discussed in 



24, 



251. When r D 



(9) 
(10) 
oo, the 



squashing function k in (2) tends to be unity so that the metric (1) reduces to that of five-dimensional RN 
black hole. The entropy and the Hawking temperature also reduce to those of five- dimensional RN cases. 
When r_ goes to zero, we have the metric for the neutral black hole with squashed horizon. If we have 
r_ = r_|_, we have deformed five-dimensional extremal RN black hole with one degenerate horizon whose 
Hawking temperature vanishes. For the case r + ,r_ — > r^, with p± finite, it is convenient to see from (6) 
that because po — > 0, K 2 (p) — ► 1, the metric describes the four-dimensional RN black hole with a twisted S 1 
bundle, where the size of the S fiber takes the constant value w p+p- = foo/2 
to that of a four-dimensional black hole in this limit. 



24j. Its temperature reduces 



The heat capacity of the black hole for the fixed Q is given by 25 



Since 3r^ — r ^ r 2 _ ^(ri + r?_) > 0, the sign of the heat capacity Cq is determined by the term r g_ (5 r 



25( and 



r +) ~ r -{^ r ^a — r +)(3 r + +?"?.) in the denominator. It was found that the Davies' point exists at r cr n 
one has Cq > for r + < r cr i t while Cq < for r + > r c rit and Cq diverges when r + crosses this critical value 
r cr i t . In the new coordinate, the Davies' point of the divergence of heat capacity Cq can be obtained from 



Pa(p+ - 5p_) + (p+ - 'ip-){p+ + p-) = 0. 

If we take b = p + — p_ , the above equation (TT2")) can be rewritten as 

ft = (P+ + P-)(P+ +P-+ 2po) 



(12) 



(13) 



2(p+ + p-) + 3 Pa 

In the limit po — * 0, we have b = (p+ + p-)/2 and p + = 3/3-, which is the Davies' point of thermal stability 
in the four-dimensional RN black hole. When pa — > oo, i.e., roc, — > oo, we obtain 6 = 2/3(p+ + /?_) and 
/?_)_ = 5p_, which is consistent with the Davies' point in five-dimensional RN black hole. 

In the following we are going to investigate whether the Davies' points on thermal stability can be reflected 
in QNM. We will concentrate on the massless scalar perturbation around the charged KK black hole with 
squashed horizons. The wave equation for the massless scalar field $(t, p, 9, <f>, ip) m the background (|6|) obeys 

1 



-9 



0p(v=5s'" , ao*(T- J p,M,v) = o. 



(14) 



Taking the ansatz $(r, p, 6, 4>, ip) = e lUT R{p)e 
we can obtain the equation 



imtfi-\-i\ip 



S(6), where S(9) is the so-called spheroidal harmonics, 



1 d 



sin0- 



dO 



S(6) 



(m — X cos ( 
sin 2 6 



E, 



lrn\ 



S(0) = o, 



(15) 



for the angular part. The eigenvalue of this angular equation (p~5|) is Ei m \ = L(L + 1) — A 2 . The radial 
equation reads 



F{P) 



p 2 K 2 {p) dp 



P 2 F{ P ) 



dR{p) 
dp 



+ [u 2 - V(p)]R(p) = 0, 



(16) 



with 



V(p) = F(p) 



L(L + 1)-X 2 4X 2 K 2 (p) 



(17) 



p 2 K 2 (p) 

The second term in the effective potential came from the fifth dimension of the spacetime. It plays a role 
of the mass of the field in the radial equation. In general, due to the presence of this term, it is difficult to 
calculate the QNM through continued fraction method. 



6 



Boundary conditions on the wave function R(p) at the outer horizon and the spatial infinity can be expressed 



as 



■ 3 / 2 / — i — 



(p-p+) "--"+ 



p ->■ p+; 



Rip) ~ { (18) 

i(P0+2(p_+p + )) ro 2 i(2p +p_+p + )A 2 
p 2 X 2(p +p_)(p +p + )x g«XP ^ p— >00. 

A solution of equation p^)l that satisfies the above boundary condition can be written as 



R(p) 



<P+ ypp + P-^ i[p +2(p_+p + )]^ 2 i( 2 P +P_+P + )A 2 



2 !po+p-Hpo+p+)x 



(19) 



(P~P+) P - _P+ XI 



m=0 



p- p+ 
p- p- 



where \ 2 = cj 2 



(po+p-)(po+p+) 

the recurrence relation starting from ao = 1 



The sequence of expansion coefficient a m : m = 1, 2, 3... is determined by 



a ai + /3 ao = 0, 

a m Om+i + P m a m + 7mffl m -i =0, m = 1, 2, 



(20) 



The recurrence coefficients a m , /3 m , 7 m are given by 



a TO = m 2 + (Co + l)m + C , 
/? m = -2m 2 + (Ci + 2)m + C 3 



7m = m 2 + (C a - 3)m + C 4 - C 2 + 2, 



(21) 



where C m are 



Co = 



2 VPo + P+^P+ 2 
P- - P+ 



P- - P+ 

«[po + 2(p- +p+)]w 2 2i / 4 /2 v /po + p+w i(2p +p- +P+)A" 

C-2 — o 1 r 7 T7 r — , 

X P--P+ IPo + P-)(Po +P+JX 



c, = 



2(po + P-)(p- - P+)(po + P+)x 
+ 2 P+ 2 Vi°o + P+(2po - P- + 3/? + )w 



(22) 



4 Z VPo + P+upT , i[( i oo + P-)(p§ + 5p +i0o + 4 i o 2 h )w 2 + (-2 i0o + /?_ -3/9+) A 2 ] 

Gl ~~ 1 1 ! v i ^ 4 ' i l6 > 

IPo + P-)(Po +P+JX 



(24) 



C 4 - 



2(po + P+)x(P- - P+) 2 + *(-2po + P- - 3/9+) (/J- - p+) 
-(Po+p-)(po + ,o+) ^2p 2 + (3po + ip+)x-i(p- -P+) (25) 

(po + 4/7+) j^ 2 + 2 / 9 + / Vpo+P+((po + 4p+)c^ 2 + 2i X ) w + 2 (i 2 + X + l) (p_ - p+)x }, 

> 3 _ | 1 f 2tp^ V/Oq + p+uj z[p + 2(/9_+/9+)]^ 2 | ^pp + p- +p+)A 2 X 
2 4\ P--P+ X (Po + P-)(Po + P+)x J 



( go +p_V 2 pg 
(P+-P-) 2 "4\- 



(26) 



If the boundary condition (1 18[) is satisfied and the series in (|20|) converge for the given L, the frequency w 
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FIG. 1: The left two panels are trajectories in the complex u) plane of scalar QNMs around squashed KK black holes, 
for n = 0, n = 1, with L = and po = 0. The others are real parts (Re(ui)) and imaginary parts (Jm(aj)) versus 
parameter a = 1 — 6. 
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FIG. 2: The left two panels are trajectories in the complex uj plane of scalar QNMs around squashed KK black holes, 
for n = 0, n = 1, with L = and po = 1- The others are real parts (Re(uj)) and imaginary parts (Im(uj)) versus 
parameter a = 1 — b. 



is a root of the continued fraction equation 



Pm—1~ Pm-2~ 



Pa 



Oimlm+l «m+l7m+2 a m+ 2"fm+3 



(m= 1,2,...). 



(27) 



Pm+l— Pm+2~ Pm+3~ 

This means that we can calculate the QNM frequencies of the charged KK black hole with squashed horizons 
by solving the above continued fraction equation (|27[). This equation is impossible to be solved analytically. 
We can only rely on the numerical calculation to obtain the QNM frequencies. 

Taking the limit Q — > 0, we have compared our numerical results by using Leaver's method with that 



obtained in 27 1 by other methods. The comparison is shown in Table HI which shows good agreement and in 
addition we see that Leaver's method we adopted gives more precise result. 

In Figs. (H}-© we display the QNM frequencies of scalar perturbation around the squashed KK black holes 



with charge for L — 0, A = 0, with overtone numbers n = 0, 1 and p — 0, 1, 2, respectively. We plot both the 
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FIG. 3: The left two panels are trajectories in the complex ui plane of scalar QNMs around squashed KK black holes, 
for n = 0, n = 1, with L = and po = 2. The others are real parts (Re(ui)) and imaginary parts (Im(ui)) versus 
parameter a = 1 — 6. 



real part and the imaginary part of QNM frequencies in functions of a = 1 — 6. Black hole horizons arc related 
to a by p + = 1 — a/2, p_ = a/2. We observe that over the critical overtone number (n c = 1 for L = in the 
scalar perturbation), both the real part and the imaginary part of QNM frequencies will behave oscillatory 
when a crosses a critical value ciq and meanwhile the complex u plan will exhibit the spiral-like shape. For 
po = 0, which is the limiting case of the four-dimensional RN black hole with a twisted S 1 bundle, we observed 
that the real part of the QNM frequency arrives at its first maximum of its oscillation approximately at the 
same position as the Davies' point = 1/2. The difference between the critical point from the QNM oq 
from that of the Davies' point is very small, |oq — ajj\ — 0.022. For po not equaling to zero, the Davies' 
thermal stability point cld can be calculated through (fT2"l) and from the behavior of QNM frequencies we can 
read the critical value oq when the oscillation of the real and imaginary parts of frequencies start. The results 
are shown in table Hi] It is interesting to find that critical points got from QNM agree very well to those of 
Davies' points. 

In Figs, (d])-©, we show results for L = 1 and A = 0, with overtone numbers n = 3, 4 and p = 0, 1, 2, 
respectively. The Davies' points are the same for each selected values of po no matter the change of L,n. As 
was observed in 17[, with the increase of L, we will observe the spiral-like shape of the complex to plan and 
oscillatory real and imaginary frequencies of QNM over higher critical overtone number n c , namely n c = 4 
when L — 1 for scalar perturbation. The oscillations of the QNM frequencies start when a over the critical 
value etQ, which is again very much in agreement with Davies' critical point ap. The results of their differences 
are shown in table [II] Choosing very big value of po, we observed that the critical point of a from QNM will 
go towards 1/3 which is the Davies point for the five-dimensional RN black hole for po — > oo. 



9 



-1.40 
-1.45 
-1.50 
-1.55 
-1.60 



(Lji)=(l,3) 






0.0 0.2 0.4 0.6 0.8 1.0 




0.38 0.40 0.42 0.44 0.46 0.48 0.50 





FIG. 4: The left two panels are trajectories in the complex lj plane of scalar QNMs around squashed KK black holes, 
for n = 3, n = 4, with L = 1 and po = 0. The others are real parts (Re(u))) and imaginary parts (Im(u))) versus 
parameter a = 1 — b. 
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FIG. 5: The left two panels are trajectories in the complex uj plane of scalar QNMs around squashed KK black holes, 
for n = 3, n = 4, with L = 1 and po = 1. The others are real parts (Re(co)) and imaginary parts (Im(uj)) versus 
parameter a = 1 — b. 
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FIG. 6: The left two panels are trajectories in the complex uj plane of scalar QNMs around squashed KK black holes, 
for n = 3, n = 4, with L = 1 and po = 2. The others are real parts (Re(uj)) and imaginary parts (Im(uj)) versus 
parameter a = 1 — b. 
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In the following we report our results for taking A ^ 0. Numerical calculation for the case A ^ is much 
more time consuming than the case with A = 0. In Fig. ((7|), we show results for L=l, n=3, 4 with po — 2, 5, 
respectively by choosing A = 0.5. In Fig. ([5]) we show results for the same choice of L, n but with A = 1. 
From these figures, it is obvious that the spiral-like shape still appears at the critical overtone number n c , 
namely n c = 4 for L = 1 with different po and nonzero A in the scalar perturbation. The critical moment to 
exhibit the spiral behavior in QNM again agrees well with the Davies thermodynamical point which is shown 
in Table HE 



TABLE I: Numerical data compared with H. Ishihara's work , with po = 0.5, p+ = 1. 



L 


A 


WKB(6th order) 


Frobrnius 


Leaver 


10 





3.5194-0.15956i 


3.510564-0.159674i 


3.5105636-0. 1596744i 


10 


1/2 


3.5336-0.15821i 


3.529225-0.158317i 


3.5292246-0. 1583174i 


10 


1 


3.5898-0. 15413i 


3.585417-0.154231i 


3.5854166-0. 1542309i 


10 


3/2 


3.6842-0. 14729i 


3.679778-0. 147368i 


3.6797783-0. 1473680i 


10 


2 


3.8178-0.13759i 


3.813421-0.137647i 


3.8134213-0. 1376469i 


10 


5/2 


3.9924-0. 12493i 


3.988012-0. 124944i 


3.98801 16-0. 1249440i 


10 


3 


4.2103-0.10912i 


4.205922-0. 109076i 


4.2059218-0. 1090755i 


10 


7/2 


4.4747-0.08988i 


4.470523-0.089745i 


4.4705233-0.0897455i 



TABLE II: The comparison of values of critical points for QNMs (a<p) and the singular point of the heat capacity (an), 
when A = 0. 



(X,L, n, po) 


Aa = \aq — an\ 


Aa 
a Q 


(0, 0, 1, 0.) 


0.022 


4.21 % 


(0, 0, 1, 1.) 


0.016 


3.85 % 


(0, 0, 1, 2.) 


0.018 


3.85 % 


(0, 1, 4, 0.) 


0.007 


0.953 % 


(0, 1, 4, 1.) 


0.0082 


2.091 % 


(0, 1, 4, 2.) 


0.0041 


1.105 % 



TABLE III: The comparison of values of critical points for QNMs (o,q) and the singular point of the heat capacity 
(od), when A = 0.5 and A = 1. 



(A, L, n, po) 


Aa = |oq — ao\ 


a O 


(0.5, 1, 4, 2.) 


0.0035 


0.942 % 


(0.5, 1, 4, 5.) 


0.0237 


6.29 % 


(1 ,1,4, 2.) 


0.0001 


0.027 % 


(1 ,1,4, 5.) 


0.0330 


8.56 % 



In summary we have studied the QNM of scalar perturbation in the background of the charged KK black 
hole with squashed horizons. We observed that over some critical value of the black hole parameter, both the 
real part and the imaginary part of the QNMs will experience oscillations and the complex to plan will exhibit 
the spiral-like shape. Interestingly this critical value agrees well to the Davies' point on the thermal stability 
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FIG. 7: the left four panels are trajectories in the complex ui plane of scalar QNMs around squashed KK black holes, 
for n = 3, n — 4, with L = 1, A = 0.5, po = 2 and 5. The others are real parts (Re(ui)) and imaginary parts (Im(ui)) 
versus parameter a. 



of black holes obtained from the singular position of the heat capacity. Our limiting case returns to the RN 
black hole, where the relation was observed in [17]. The correspondence of the critical point observed in QNM 
to the position of infinite discontinuities of the heat capacities indicates that QNM may shed the light on the 
turning point of the thermal stability. It is of interest to generalize this discussion to other backgrounds and 
for different fields' perturbations. 

For a black brane solution it was conjectured that black holes which lack local thermodynamical stability 
often also lack stability against small perturbations 23[. This conjecture might only hold for black holes with 



translation symmetry, such as black string. In our background spacetime we observed that on both sides 
of the turning point of the thermodynamical stability, the imaginary frequencies of the QNM are negative, 
which tells us that the scalar perturbation is always stable even at the critical point. Recently the metric 



perturbations were discussed 
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271 ] and it was indicated that dynamically the squashed KK black hole is 
stable, this is consistent with our result. This result shows that in the background we are studying it has 
dynamical stability even when thermodynamical instability appears. This is quite interesting, which suggests 
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FIG. 8: the left four panels are trajectories in the complex ui plane of scalar QNMs around squashed KK black holes, 
for n — 3, n — 4, with L = 1, A = 1, po = 2 and 5. The others are real parts (Re(ui)) and imaginary parts (7m(oj)) 
versus parameter a. 

that the relation between the dynamical and thermodynamical stabilities of black holes is non-trivial and 
further investigations are called for. 
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